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Abstract 

Transport phenomena are ubiquitous in nature and known to be impor- 
tant for various scientific domains. Examples can be found in physics, elec- 
trochemistry, heterogeneous catalysis, physiology, etc. To obtain new infor- 
mation about diffusive or Laplacian transport towards a semi-permeable or 
resistive interface, one can study the random trajectories of diffusing parti- 
cles modeled, in a first approximation, by the partially reflected Brownian 
motion. This stochastic process turns out to be a convenient mathematical 
foundation for discrete, semi-continuous and continuous theoretical descrip- 
tions of diffusive transport. 

This paper presents an overview of these topics with a special empha- 
sis on the close relation between stochastic processes with partial reflections 
and Laplacian transport phenomena. We give selected examples of these 
phenomena followed by a brief introduction to the partially reflected Brow- 
nian motion and related probabilistic topics (e.g., local time process and 
spread harmonic measure). A particular attention is paid to the use of the 
Dirichlet-to-Neumann operator. Some practical consequences and further 
perspectives are discussed. 
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Introduction 

An erratic motion of pollens of Clarkia (primrose family) discovered by Robert 
Brown in 1827 and quantitatively described by Albert Einstein in 1905 gave a 
substantial impact for developing mathematical theory of stochastic processes, an 
important branch of modern mathematics. Supported by rigorous mathematical 
foundations, the Brownian motion and related stochastic processes found numer- 
ous applications in different scientific domains, from theoretical physics to biol- 
ogy and economics. To study the transport of species diffusing from a remote 
source towards and across semi-permeable or resistive interface (e.g., cellular 
membrane), one can employ either an averaged description in terms of an ap- 
propriate boundary value problem for the concentration of species, or stochastic 
analysis of their random trajectories. In the first case, a finite permeability (re- 
activity, resistivity, etc.) of the interface leads to the mixed or Fourier boundary 
condition, while in the second case it can be modeled as partial reflections on the 
boundary. Physical or chemical processes governed by the Laplace equation (sta- 
tionary diffusion) with mixed boundary condition are generally called Laplacian 
transport phenomena. Their examples are found in physiology (oxygen diffusion 
towards and across alveolar membranes), in electrochemistry (electric transport 
in electrolytic cells), in heterogeneous catalysis (diffusion of reactive molecules 
towards catalytic surfaces), in nuclear magnetic resonance (diffusion of spins in 
confining porous media). Studying random trajectories of diffusing species, one 
can extract a subtle information about the system in question. However, the theo- 
retical or numerical analysis of these phenomena is in general complicated by an 
irregular geometry of the interface (e.g., microroughness of metallic electrodes, 
see Section [T]). 

In this paper, we focus on a particular stochastic process, called partially re- 
flected Brownian motion (PRBM), and its application to study Laplacian transport 
phenomena. Our main purpose is to capture the attention to this interesting pro- 
cess itself, and its use for understanding the influence of a geometrical irregularity 
of the interface on transport properties. Bearing in mind the particular role of the 
geometry, we would like to "bridge" theoretical, numerical and experimental stud- 
ies of Laplacian transport phenomena on the one side, and powerful mathematical 
methods of stochastic analysis on the other side. Since the extensive literature 
existing on both topics is generally difficult to get through for non-specialists, we 
prefer to use a descriptive style of writing in order to give the whole vision of the 
problem, without specifying particular details which can be found anywhere else 
(e.g., see references at the end of this paper). 

In the first section, we present three examples of Laplacian transport phenom- 
ena in different scientific fields. Their mathematical description by the mixed 
boundary value problem opens encouraging possibilities to apply powerful tools 



Partially Reflected Brownian Motion. 



3 



of potential theory, variational analysis and probability theory. The second section 
is devoted to remind some basic definitions of stochastic process theory: stop- 
ping times, reflected Brownian motion, local time process, harmonic measure, 
etc. In the third section, we introduce the partially reflected Brownian motion and 
show its properties for a planar surface. An important relation to the Dirichlet- 
to-Neumann operator is revealed and then illustrated by several examples. The 
last section presents different stochastic descriptions of Laplacian transport phe- 
nomena: a recently developed continuous approach and two other methods. In the 
conclusion, we summarize the essential issues of the paper. 

1 Laplacian Transport Phenomena 

The transport of species between two distinct "regions" separated by an interface 
occurs in various biological systems: water and minerals are pumped by roots 
from the earth, ions and biological species penetrate through cellular membranes, 
oxygen molecules diffuse towards and pass through alveolar ducts, and so on. 
Transport processes are relevant for many other scientific domains, for example, 
heterogeneous catalysis and electrochemistry. In this section, we shall give three 1 
important examples of the particular transport process, called Laplacian or diffu- 
sive transport. 

1.1 Stationary Diffusion across Semi-permeable Membranes 

Let us begin by considering the respiration process of mammals. Inbreathing a 
fresh air, one makes it flow from the mouth to the dichotomic bronchial tree of the 
lungs (Fig.[l]). For humans, first fifteen generations of this tree serve for convec- 
tional transport of the air towards pulmonary acini, terminal gas exchange units 
JH]. A gradual increase of the total cross-section of branchiae leads to a decrease 
of air velocity. At the entrance of the acinus, it becomes lower than the character- 
istic diffusion velocity ||2|]. As a consequence, one can describe the gas exchange 
inside the acinus as stationary diffusion of oxygen molecules in air from the en- 
trance ("source" with constant concentration Co during one cycle of respiration) 
to the alveolar membranes 0]. In the bulk, the flux density is proportional to the 
gradient of concentration (Fick's law), J = —DVC, where D is the diffusion 
coefficient. The mass conservation law, written locally as div J = 0, leads to the 
Laplace equation AC = in the bulk. The flux density towards the interface is 
simply J n = D dC/dn, where the normal derivative d/dnis directed to the bulk. 

'Diffusive NMR phenomena present another important example when the transport properties 
are considerably affected by irregular geometry. In this paper, we do not discuss this case since it 
has been recently reviewed in a separate paper [4]. 
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Arrived to the alveolar membrane, oxygen molecules can penetrate across the 
boundary for further absorption in blood, or to be "bounced" on it and to continue 
the motion. The "proportion" of absorbed and reflected molecules can be charac- 
terized by permeability W varying from (perfectly reflecting boundary) to in- 
finity (perfectly absorbing boundary). In this description, the flux density across 
the alveolar membrane is proportional to the concentration, J n = WC. Equat- 
ing these two densities on the alveolar membrane, one gets a mixed boundary 
condition, D(dC/dn) = WC, called also Fourier or Robin boundary condition. 
Resuming these relations, one provides the following mathematical description 
for the diffusion regime of human or, in general, mammalian respiration: 



A— 

dn 



AC = in the bulk (1) 
C = Co on the source (2) 

C = on the alveolar membrane (3) 



where the underlying physics and physiology are characterized by a single pa- 
rameter A = D/W, which is homogeneous to a length (/ stands for the identity 
operator). Note also that the dependence on constant C is irrelevant. In what 
follows, we address to this "classical" boundary value problem. The essential 
complication resides in a very irregular geometry of the pulmonary acinus, which 
presents a branched structure of eight generations (for humans), "sticked" by alve- 
olar ducts (Fig.[T]). For small A, only a minor part of the boundary is involved to 
the transport process (so-called Dirichlet active zone), whereas the flux across 
the rest of the boundary is almost zero (this effect is called diffusional screening 
13, El EELS, El)- With an increase of A, larger and larger part of the boundary be- 
comes active. As a result, the efficiency of human lungs depends on the parameter 
A in a nontrivial manner that implies different physiological consequences Jlflll . 
The trajectory of a chosen oxygen molecule can be seen as Brownian motion from 
the source towards the alveolar membrane, with multiple bounces on the boundary 
and final absorption. This is in fact what we call the partially reflected Brownian 
motion (Section A profound study of the interplay between the irregular ge- 
ometry of the acinus and the erratic random motion of oxygen molecules inside it 
should help to better understand physiological functioning of human lungs. 



1.2 Heterogeneous Catalysis 

A similar description can be brought to the molecular regime of heterogeneous 
catalysis omnipresent in petrochemistry. One considers reactive molecules A in- 
jected into a solvent and then diffusing towards a catalyst. Hitting the catalytic 
surface, they can be transformed into other molecules A* (with a finite reaction 
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Figure 1: On the left, a cast of human lungs; on the right, a cast of pulmonary 
acinus (HQ (by E. Weibel). 



rate K), or to be bounced for further diffusion in the bulk. The new molecules 
A*, collected by appropriate physical or chemical technique, do not further con- 
tribute to the transport process. Assuming the presence of a remote source of 
reactive molecules A, one can model, in a first approximation 2 , the heterogeneous 
catalysis by the mixed boundary value problem ((l]-|3l) with a characteristic length 
A = D/K JiEEEE3,El]. The keynote of this similitude is related to the fact 
that each reactive molecule arrived onto the boundary terminates its motion af- 
ter a number of successive reflections. The mechanism leading to its termination 
is different: for the oxygen diffusion, the molecules are absorbed by the alveo- 
lar membrane and transferred to the blood, while for the heterogeneous catalysis, 
the reactive molecules are transformed by chemical reaction into other molecules 
which do not further participate to the process. Since the overall production of 
new molecules A* depends on the total surface area of the catalytic surface, one 
tries to design catalysts with the largest possible surface (for given volume), real- 
izing porous and very irregular boundaries (Fig. El). As a consequence, the diffu- 
sional screening becomes important to understand numerous industrial processes 
in petrochemistry. Since random trajectories of reactive molecules correspond to 
the partially reflected Brownian motion, its study may allow a design of more 



2 This description is probably too simplified in order to model the heterogeneous catalysis quan- 
titatively. First, the presence of molecules A* near the catalyst may "obstruct" the access to the 
catalytic surface. Second, parasite reactions happen on the boundary that implies a progressive 
deactivation of the catalyst. Consequently, the reactivity K becomes dependent on the spatial 
position on the catalytic surface, leading to an inhomogeneous boundary condition. Finally, the 
molecular diffusion can be applied only if the mean free path of reactive molecules is much lower 
than the geometrical features of the catalyst (in the opposite case, one deals with Knudsen dif- 
fusion 112, III El LLJ). Nevertheless, the simple description QH) permits to take into account 
many important features related to the catalytic process. 
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Figure 2: On the left, an example of an irregular catalytic surface (by J. S. An- 
drade jr.); at the center, photo of a rough metallic surface of nickel electrode (by 
E. Chassaing); on the right, photo of an irregular metallic electrode used to study 
the Laplacian transport phenomena experimentally (by B. Sapoval). 



efficient catalysts. 



1.3 Electric Transport in Electrochemistry 

The other example of Laplacian transport phenomena can be found in electro- 
chemistry: the electric current between two metallic electrodes into an electrolyte 
is described by the same boundary value problem. Indeed, the electric potential V 
obeys the Laplace equation in the bulk since the electrolyte is locally neutral. Tak- 
ing one electrode of very low resistance (counter-electrode), one writes the corre- 
sponding boundary condition as V = V , where V is the applied tension. For the 
other electrode of surface resistance r (working electrode), one obtains the mixed 
boundary condition by equating the volume current density — p~ l W (p is the 
electrolyte resistivity) and the surface current density V/r: A dV/dn = V, where 
A = r/p is again the physical length of the problem. The similar description can 
be brought even in the case of an alternative tension [ 19, 20]. 

For electric transport, one cannot associate directly the mixed boundary value 
problem with the partially reflected Brownian motion since there is no diffus- 
ing particle. From this point of view, the electrochemical problem has only a 
formal analogy with two previous examples. At the same time, the electrochem- 
istry is an appropriate domain to study experimentally the influence of the ir- 
regular geometry on the (average) transport properties. Taking metallic elec- 
trodes of different shapes with micro- or macro-roughness (e.g., see Fig.EJ), one 
can directly measure the spectroscopic impedance or admittance (see below). 
These characteristics are equivalent to the total flux across the boundary for dif- 



fusional problems 112 lL 12211 . The observation of anomalous impedance behav- 



ior 1 23] had provoked numerous theoretical, numerical and experimental stud- 
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ies of the role of a geometrical irregularity in Laplacian transport phenomena 
H H HI H BSHH0SHHBBH(for more informa- 
tion, see [39] and references therein). 



1.4 Discrete and Semi-continuous Approaches 

Among different theoretical approaches developed to study Laplacian transport 
phenomena, we have to mention the double layer theory of Halsey and Leibig 



f— i 



[26, 27, 28, 29] and the formalism of the Brownian self-transport operator pro- 
posed by Filoche and Sapoval I2U1 . In Section l4~3l we shall show how the origi- 
nal Green function description by Halsey and Leibig can be related to the Brow- 
nian motion reflected with jump (this stochastic reformulation will be referred to 
as "semi-continuous" approach). In turn, Filoche and Sapoval considered lattice 
random walks with partial reflections to derive a spectral representation for the 
macroscopic response of an irregular interface (see Section 14.41 where this for- 
malism is referred to as "discrete" approach). Although both methods accurately 
describe Laplacian transport phenomena (e.g., they give an explicit formula for 
the total flux across the boundary), their major inconvenience resides in the de- 
pendence on an artificial length scale: jump distance a for the semi-continuous 
approach and lattice parameter a for the discrete approach. A physical intuition 
suggests that, if these descriptions are correct, there should exist a well defined 
continuous limit as a tends to 0. Certain substantial arguments to justify the ex- 
istence of this limit were brought in IB9I1 (and they will be strengthened in this 
paper), but a rigorous mathematical proof is still required. To overcome this dif- 
ficulty, a new theoretical approach has been recently developed in lEoll . We shall 
call it "continuous" since it is tightly related to a continuous stochastic process, 
namely, the partially reflected Brownian motion. This approach will integrate 
the advantages of the previous ones, being a mathematical foundation for under- 
standing Laplacian transport phenomena. We shall return to these questions in 
Section El 



2 Basic Definitions 



In this section, we recall the basic definitions related to the Brownian motion and 
reflected Brownian motion that can be found in extensive literature, e.g., H41LI42 . 
El . 44 , 4J , 4^, 42 ] . The familiar reader may pass over this section. 
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2.1 Brownian Motion and Dirichlet Boundary Value Problem 



The Brownian motion can be defined in different ways 114 111 . Throughout this 
paper, we use the following definition. 

Definition 2.1 A stochastic process W t (t > 0) defined on the chosen probabilis- 
tic space is called one-dimensional Brownian motion (or Wiener process) started 
from the origin, if 

• its trajectories are continuous almost surely (with probability 1); 

• it starts from the origin almost surely, P{Wo = 0} = 1; 

• its joint distribution is 

F{W tl eT 1 ,...,W tn eT n } = 

tfei... / dx n g(0, x 1 ; h) g(x u x 2 ; t 2 - h) ... g{x n ^, x n ; t n - t n -i] 



for any integer n, any real numbers < t\ < ... < t n and arbitrary inter- 
vals T lt r n , where g(x, x' ; t) is the Gaussian density 



g(x,x' ; t) 



'2id 



cxp 



[X — X 



2t 



x, x e R, te M+ (4) 



By definition, g(0, x ; t)dx is the probability to find the Brownian motion in 
dx vicinity of point x at time t: 

F{W t e (x, x + dx)} = g(0, x ; t)dx 

The collection W t = (W/, Wf) of d independent one-dimensional Brown- 
ian motions is called d-dimensional Brownian motion started from the origin 
( in the following, we shall omit the pointing on the dimension ). The translated 
stochastic process, x + W t , is called Brownian motion started from the point 
x E R d . 

Various properties of the Brownian motion and its relation to other mathemat- 
ical fields (like partially differential equations or potential theory) are well known 
and can be found in (Jj 51 BAB EH. 

As one can see, the Brownian motion W t is defined for the whole space R. d , 
without any binding to a particular domain. However, physical processes are usu- 
ally confined into a certain domain f2 C W 1 . The "presence" of its boundary 
<9f2 can be introduced by a specific condition for a quantity we are looking for. 
To illustrate this notion, let us introduce the harmonic measure cu x defined as the 
probability measure to hit different subsets of the boundary dtt for the first time. 
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Definition 2.2 Let Q C M> d be a domain with boundary dQ. For any x G M, d , a 
random variable TP = inf{t > : (x + W t ) G dVt] is called stopping time 
on the boundary dVt (it gives the first moment when the Brownian motion started 
from x hits the boundary). For any subset A from the Borel a -algebra B(dQ), one 
defines its harmonic measure lo x {A} (hitting probability) as: 

u x {A} = F{W T * G A, TP < oo} 

(we remind that the Borel a -algebra B(dQ) is generated by all open subsets of 
dQ). 

We gave this classical definition of the harmonic measure in order to outline 
that the boundary dVl is present in the problem only through the stopping time TP. 
In other words, its introduction does not change the definition of the Brownian 
motion itself. This feature considerably simplifies the following analysis. 

Up to this moment, we did not specify the domain and its boundary d£L 
since the harmonic measure can be well defined for very irregular domains [48, 
I49I lioll . However, the following definitions will need some restrictions on the 
boundary. Throughout this paper, we shall consider a domain C W 1 (d > 2) 
with bounded smooth boundary dVl (twice continuous differentiable manifold). 
One the one hand, this condition can be weakened in different ways, but it would 
require more sophisticated analysis overflowing the frames of this paper (e.g., see 
El El). One the other hand, our primary aim is to describe Laplacian transport 
phenomena listened in Section U Dealing with physical problems, one can al- 
ways smooth a given boundary d£l whatever its original irregularity. Indeed, the 
physics naturally provides a minimal cut-off 5 (e.g., mean free path of diffusing 
or reacting molecules) which determines the "admissible" scales of the bound- 
ary. All geometrical features of the boundary smaller than 5 should be irrelevant 
(otherwise, the proposed physical description would be incorrect). Smoothing 
these geometrical elements, one can obtain a boundary that may be (very) irregu- 
lar on length scales larger than 5, but smooth on length scales lower than 5. For 
a smooth boundary d£l, one can introduce the harmonic measure density u> x (s) 
such that co x (s)ds is the probability that the Brownian motion started at x hits the 
boundary in ds vicinity of the boundary point s. 

The harmonic measure, generated by the Brownian motion, gives a general 
solution of the Dirichlet boundary value problem with a given function / on dVt: 

Au = (x g n), u = f (xe on) (5) 

Indeed, the harmonic measure density is equal to the normal derivative of the 
Green function for the Dirichlet problem, so that one writes the solution u(x) 
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explicitly Q: 

u(x) = J f(s) u x (s)ds 
an 

or as following expectation 

u(x) = E{f(W T *)} (6) 

One can give a physical interpretation to this mathematical relation. In order to 
calculate the expectation, one considers all possible trajectories of the Brownian 
motion started from the point x G fi. For each trajectory terminated at boundary 
points s = Wjx, one assigns the weight f(s) and then averages over all these 
trajectories. Giving this interpretation, we do not discuss the mathematical real- 
ization of such average over all possible trajectories. To do this operation properly, 
one can introduce the Wiener measure on the space of continuous functions and 
then define the corresponding functional integrals 114 ill . It is interesting to remark 
that this reasoning traced to the Feynman's description of quantum mechanics by 
path integrals 115411 . Note also that the relation © is the mathematical foundation 
to Monte Carlo numerical tools for solving the Dirichlet problem ©. In fact, 
launching a large number of random walkers from the point x, one determines, 
for each trajectory, its hitting point s and assigns the corresponding weight f(s). 
The average over all random walkers gives an approximate value of the solution 
u(x) at point x. 

One can see that the Brownian motion is an efficient mathematical tool to 
study Dirichlet boundary value problems. However, it becomes useless for other 
types of boundary conditions like, e.g., the Neumann condition. The simple phys- 
ical reason is the following. As we have mentioned above, the Dirichlet boundary 
condition is introduced through the stopping time T x . It means that we are inter- 
ested only in the Brownian motion W t for times t between and T x . Since the 
motion with t > T x is irrelevant for this problem, one may think that the Brownian 
motion is absorbed on the boundary d£l at the first hit. In other words, the Dirich- 
let condition corresponds to a purely absorbing interface d£l. For the Neumann 
condition, the situation changes drastically. The normal derivative representing 
a flux leads to the notion of reflection on the boundary: if one would like to fix 
the flux density across the boundary, certain particles should be reflected. The 
probabilistic description of the Neumann boundary condition necessitates thus an 
introduction of the other stochastic process called reflected Brownian motion. 

2.2 Reflected Brownian Motion 

The fact of reflection on the boundary implies three essential distinctions with 
respect to the (simple) Brownian motion: 
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• the definition of the reflected Brownian motion will depend on the domain f2 
(as a consequence, it will be necessarily more sophisticated than the above 
definition of the Brownian motion); 

• the type and direction of each reflection should be prescribed (e.g., right or 
oblique); 

• some restrictions on the boundary dVl should be introduced, for instance, 
the normal vector should be well defined at each point (as a consequence, 
the boundary cannot be very irregular). 

It is not thus surprising that the definition of the reflected Brownian motion re- 
quires stochastic differential equations. We do not intend to reproduce the whole 
analysis leading to the reflected Brownian motion since one can find it in corre- 
sponding literature (e.g., see J4lll55 . 56i 5t\\ ). In the case of smooth boundaries, 



the following definition is quite classical. The situation becomes essentially more 
difficult when one tries to extend it for nonsmooth domains. 

Definition 2.3 Let Q C M. d be a domain with boundary dQ, and n(s) is a vector- 
valued function on dfl. For a given point x G Q, one considers the stochastic 
equation in the following form: 

dW t = dW t + n(W t )I dn (W t )d£ t W = x, £ = (7) 

where W t is d-dimensional Brownian motion and Iqq is the indicator of the bound- 
ary dfl. By a solution of this equation, we mean a pair of almost surely continuous 
processes W t and £ t , satisfying (0), adapted to the underlying family of a -fields 
and satisfying, with probability 1, the following conditions: 

• W t belongs toQU dfl; 

• £ t is a nondecreasing process which increases only for t G T, T = {t > 
: W t G dQ} having Lebesgue measure zero almost surely. 

The process W t is called Brownian motion normally reflected on the boundary ( or 
reflected Brownian motion), the process it is called local time on the boundary 
(or local time process). 

The following theorem ensures the existence and uniqueness of these stochas- 
tic processes in the case of smooth boundaries. 

Theorem 2.4 Let Q C M d be a bounded domain with twice continuous differ- 
entiable boundary dQ, n(s) is the vector of the inward unit normal at boundary 
point s ( orthogonal to the boundary at s and oriented towards the domain ). For 
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a given point x G Q U dQ, the stochastic equation (0) possesses a unique solu- 
tion, i.e., there exist the reflected Brownian motion Wt and the local time on the 
boundary £ t satisfying the above conditions, and they are unique. 



Proof can be found in 14 lL 15511 



We should note that this theorem can be extended in different ways. For ex- 
ample, one can consider the Brownian motion, reflected on the boundary in the 
direction given by another vector- valued field than the field n(s) of the inward 
unit normals. The assumption that the domain is bounded can be replaced by a 
more subtle hypothesis that allows to extend the definition of the reflected Brow- 
nian motion for some classes of unbounded domains. At last, one may define 
this motion for a general case of second order elliptic differential operators (with 
certain restrictions on their coefficients). The interested reader may consult the 



corresponding literature, e.g., [4fL 14511 . 

Although the rigorous mathematical definition of stochastic differential equa- 
tions is more difficult than in the case of ordinary differential equations, an in- 
tuitive meaning of its elements remains qualitatively the same. For example, the 
stochastic equation © states that an infinitesimal variation dW t of the reflected 
Brownian motion W t in the domain (bulk) is governed only by the variation 
dW t of the (simple) Brownian motion W t (the second term vanishes due to the in- 
dicator Ida). When the motion hits the boundary, the second term does not allow 
to leave the domain leading to a variation directed along the inward unit normal 
n(s) towards the interior of the domain. On the other hand, each hit of the bound- 
ary increases the local time £ t . Consequently, the single stochastic equation © 
defines simultaneously two random processes, W t and £ t , strongly dependent each 
of other. 

As an example, one can consider one-dimensional Brownian motion reflected 
at zero which can be written as mirror reflection of the (simple) Brownian motion: 
W t — \x + W t \. Applying Ito's formula to this function, one obtains: 

t t 
W t = \x + W t \=x + J sign(x + W v )dW v + \ J 5{x + W t <)dt' 



One can show that the second term is equivalent to a Brownian motion W[, whereas 
the third term, denoted as £ t , is a continuous, nondecreasing random process which 
increases only on the set T = {t > : x + W t = 0} of the Lebesgue measure 
zero. The previous expression can thus be written as Wt — x + W t + £ t or, in 
differential form, as dW t = dW t + d£ t which is the particular case of the stochas- 
tic equation (JT]). For the local time £ t , Levy proved the following representation 
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EH: 

t 

£ t = UmJ- f I [0M (W t >)dt' (8) 



This relation makes explicit the meaning of the local time £ t : it shows how "many 
times" the reflected Brownian motion passed in an infinitesimal vicinity of zero 
up to the moment t. Levy also gave another useful representation for the local 
time: 

£ t = lim aXt(a) (9) 

a— >0 

where A/f (a) is the number of passages of the reflected Brownian motion through 
the interval [0, a] up to the moment t. If one introduces a sequence of stopping 
times at points and a, 

r (0) = inf {* > : W t = 0} r (a) = inf {t > r (a) : W t = a} 

rf = inf {t > : W t = 0} r« = inf {t > : W t = a} 

the number of passages can be defined as 

M t {a) = sup{n > : rf } < t} 

Note that the representations dU) and © can be extended for a general case of 
(i-dimensional reflected Brownian motion. 



3 Partially Reflected Brownian Motion 
3.1 Definition and Certain Properties 

Bearing in mind the description of Laplacian transport phenomena, we would like 
to extend the concept of the reflected Brownian motion in order to deal with the 
mixed boundary condition ©. 

Definition 3.1 For a given domain Q C M. d with smooth bounded boundary dQ, 
let Wt be the reflected Brownian motion started from x G £1 U dQ, and it be the 
related local time process. Let x be a random variable, independent of W t and £ t 
and distributed according to the exponential law with a positive parameter A: 

P{ x > A} = exp[-A/A] (A > 0) (10) 

The stopping time 

Tl = inf{t > : £ t > X} 

gives the first moment when the local time process £ t exceeds the random variable 
X- The process W t conditioned to stop at random moment t = T\ is called 
partially reflected Brownian motion (PRBM). 
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First of all, we stress that the partially reflected Brownian motion is not a new 
stochastic process: it reproduces completely the reflected Brownian motion W t up 
to the moment T^. The only difference between them resides in the fact that we 
are not interested in what happens after this moment. Consequently, the condition 
to stop dXt = T\ may be thought as an absorption on the boundary dVt. It ex- 
plains the term "partially reflected": after multiple reflections, the process will be 
absorbed on the boundary (see Section|4]for further comments). Roughly speak- 
ing, the whole term "partially reflected Brownian motion" is a shorter version of 
the phrase "reflected Brownian motion conditioned to stop at random moment 

In the particular case A = 0, the exponential distribution (TTOb is degenerated: 
^{x — 0} = 1 and F{x > 0} = 0. Consequently, the stopping time becomes: 
Tg = infjt > : i t > 0}. Since the first moment of an increase of the local 
time process £ t corresponds to the first hit of the boundary dVl, one obtains the 
stopping time of the (simple) Brownian motion: Tq = T x . One concludes that, 
for A = 0, the partially reflected Brownian motion becomes the Brownian motion 
conditioned to stop at the first hit of the boundary. 

To study the partially reflected Brownian motion, one can introduce a measure 
quantifying absorptions on different subsets of the boundary <9fl 

Definition 3.2 For any subset A from the Borel a -algebra B(dQ), one defines its 
spread harmonic measure u Xj &{A} as: 

u xA {A} = F{W n eA,T%< oo} 

As the harmonic measure itself, uj x> k{A] satisfies the properties of a prob- 
abilistic measure, in particular, uj x ^{d£i} = 1. When A goes to 0, the spread 
harmonic measure tends to the harmonic measure: co X: \{A} — ► ^{A}. Since 
the present definition of the PRBM requires the smoothness of the boundary, the 
spread harmonic measure can be characterized by its density u x j^(s). 

Dealing with the Brownian motion, one could formally take the starting point 
x on the boundary <9f2, but it would lead to trivial results: the stopping time T x 
becomes and the harmonic measure lo x is degenerated to the Dirac point mea- 
sure: c<j x .{A} = Ia(x) (if x G dil). In the case of the partially reflected Brownian 
motion, the starting point x can belong to the domain f2 or to its boundary d£l: in 
both cases the spread harmonic measure has nontrivial properties. 

It is convenient to separate each random trajectory of the PRBM in two parts, 
before and after the first hit of the boundary. The first part, Wo<t<r°>, coincides 
with the (simple) Brownian motion started from x and conditioned to stop on the 
boundary, while the second part, W^<t<T x , coincides with the reflected Brow- 
nian motion started on the boundary (at the first hitting point) and conditioned 
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to stop on the same boundary at random moment T A . Since these two parts are 
independent, one can write the spread harmonic measure density as 



w x ,a(s) = / ds'uj x (s') Ta(s' , s) 



T A (s', s) = uv,a(s) 



(11) 



i)0 



The integral kernel T\(s',s) represents the probability density that the PRBM 
started from the boundary point s' is stopped (absorbed) in an infinitesimal vicinity 
of the boundary point s. Consequently, it is sufficient to determine the probabil- 
ities of displacements between two boundary points in order to reconstruct the 
whole spread harmonic measure density. 

Lemma 3.3 For any subset A from B(dQ) and fixed positive A, the spread har- 
monic measure ui x ^{A}, considered as a function ofx, solves the mixed boundary 
value problem: 



Acu xA {A} = (xeQ), 



dn 



u xA {A} = l A {x) (xedn) 

(12) 



This lemma generalizes the Kakutani theorem for the harmonic measure (when 
A = 0) S. We do not reproduce the proof of this lemma since it would require 
many technical details. It can be also reformulated for the spread harmonic mea- 
sure density: 

Lemma 3.4 For any boundary point s e dQ and fixed positive A, the spread 
harmonic measure density uj x a(s), considered as a function of x, satisfies the 
following conditions: 



Au xA (s) = (x e Q), 



I- A 



d 



On 



(13) 

where 5(s — x) is the Dirac function (distribution) on the boundary. 



According to this lemma, the solution of a general mixed boundary value prob- 



lem 



Au = (xett), 



I-A 



d_ 

dn 



u = f (xedQ) 



with a given function / on d£l and fixed positive A can be written in two equivalent 
forms: 



u(x) 



- J f(s) u xA (s)ds = E{f(W ri )} 
an 
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Again, one can give a physical interpretation of this relation: one averages the 
function / over all possible trajectories of the partially reflected Brownian motion 
started from the point x. Each trajectory is weighted by f(s) according to the 
boundary point s of its final absorption. 

3.2 Planar Surface 

We remind that the physical motivation of this work is a possibility to describe dif- 
fusing particles near semi-permeable interfaces by the partially reflected Brownian 
motion. Indeed, the mixed boundary value problem (JTJ-El) is an averaged descrip- 
tion for the concentration of particles, while the stochastic description permits to 
"follow" the trajectory of one individual particle. This analysis may provide a 
new information: typical or average distance between the first hitting point and 
the final absorption point; proportion of "flatten" trajectories, going near the inter- 
face, with respect to remote trajectories, moving away from the interface and then 
returning to it, etc. In this subsection, we briefly consider the particular case of the 
planar surface (boundary of a half space), when the partially reflected Brownian 
motion can be constructed in a simple way, without stochastic equations. Conse- 
quently, many related characteristics can be obtained explicitly. In addition, this 
construction for the half space brings an example of the PRBM for an unbounded 
domain. 

Let f2 be the upper half space, = {x E M d : Xd > 0}, with smooth 
boundary dVl = {x E K d : Xd = 0}. Let {W 7 ^} are d independent Brownian 
motions started from the origin. Then, the Brownian motion, started from a given 
point x E fl, and reflected on the boundary <9f2, can be written in a simple way as 
(x± + Wf, Xd-i + W[ , \xd + Wf\). The particular simplification is brought 
by the fact that reflections happen in a single direction, being involved through the 
one-dimensional reflected Brownian motion \x^+ Wf\. Without loss of generality, 
we can consider the reflected Brownian motion started from the origin (x = 0): 
the translational invariance along the hyperplane d£l permits to move the starting 
point in dVl, whereas the convolution property (fTTTl allows displacements in or- 
thogonal direction. The local time process i t can be introduced either through the 
stochastic equations ©, or with the help of Levy's formulae ® or ©. 

Lemma 3.5 Let Q be the upper half space, Q = {x E M. d : Xd > 0}. For any 
positive A, the stopping time T^, defined in \3.1\ is distributed according to 




CO 



P{T' 




(14) 



o 
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Proof . Since the local time £ t and the random variable x are independent, one can 
write the probability P{T A e(t,t + dt)} as 

oo 

P{ T° A G (t, t+dt) } = J P j inf {r > : t T = z) G (t, t+dt) J P{ X G (-2, z+dz) } 



P^ inf {r > : £ T = z) G (t, t + dt) 



dt 



Then, the first factor is the well known density of the inverse local time process 

z e~ z2 l 2t 

while the second factor is given by the exponential law density (ITOb that implies 

<m. □ 

Note that the integral in (THl) can be represented with the help of the Gaussian 
error function 



PA(t) 



2A 2 



1 1 

1 v^/W 



/q vty2A 2 



One finds the asymptotic behavior of the density pa(£): 

p A (t) ~ (v^Af 1 r 1 ' 2 (t - 0), p A (t) ~ (v^/A)" 1 1- 3 / 2 (t - oo) 

Once the distribution of stopping time T A is determined, one can calculate the 
spread harmonic measure density lo X) a{s). 

Lemma 3.6 Let Q be the upper half space, Q = {x G M. d : Xd > 0}. For any 
positive A, the spread harmonic measure density cu x> a(s) is 



oo oo 



dk\...dkd- 
72 



ex P 



-oo — oo 



d-1 



-a; d |fc| 



1 + A|Jfe| 

(15) 



where \k\ 



k\ + ... + h\_ v 



Proof . First, the probability kernel T\(s,s'), defined for two boundary points 
s, s' G <9f2, is translationally invariant in the hyperplane <9f2, T\(s, s') = t\(s—s'), 
where 



tA(sx,—,Sd-i)ds 1 ...ds d -i = 

p| W^o a G (s^sx + dsx) 



, ... , H^o 1 G (s d -i, s<i-i + dsd-i) 
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The stopping time T° A is related to the orthogonal motion and, consequently, inde- 
pendent of lateral motions W}, Wf~ l . Therefore, the above probability can be 
written as 



t A (si, Sd-x)dsx...dsd-i 



oo 

j P j W] G ( Sl , Si + dSi), ... , Wf- 1 G (S,,-!, + d Sd ^) | 



P\(t)dt 



Since the lateral motions are independent between themselves, the first factor is 
equal to the product of Gaussian densities ©: 



t A (si, ...,s d _i) = j dtp A (t) Yl 
o 

Using the integral representation (TBI) , one finds 



- , v 2?rt 



, r(d/2) /• , * e-*/ A 

■■■■ «-0 = ^ ./ * H + ... + s ;_ l + ^ 







where T(z) stands for Euler gamma function (see lliill for details). 

Substituting the well known harmonic measure density lo x (s) for the upper 
half space (generalized Cauchy distribution), 

, v r(d/2) 



^ [(Xi - Sl ) 2 + ... + - S ,_0 2 + (X d ) 2 ]" /2 



into convolution (fTTT) . one finally obtains the expression (IT31) for the spread har- 
monic measure density. □ 

One can easily verify that the spread harmonic measure density u X! \( s ) an d 
the probability kernel T A (s, s') satisfy the following conditions: 

1 . Normalization condition: 

J ds u X;A {s) = 1 J ds' T A (s, s') = 1 

an an 

2. Dirichlet limit (A -> 0): 

W«b,a(s) ► Was(s) T A (s, s') >5(S- s') 
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3. Translational invariance: 

u x ,a( s ) = ^-s,a(°) t a( s > s ') = Ta(s - s', 0) = t\(s - s') 

One can also deduce the asymptotic behavior of the function t^s) as \s\ — > or 
\s\ — > oo. For this purpose, it is convenient to define the new function r]d(z) by 
relation 

t\(s) = Vd{\s\/A) W( ,... ) o,A)(s) 

where the second factor is the harmonic measure density for the Brownian mo- 
tion started from the point (0, 0, A). Using the explicit formulae for t\(s) and 

d-l 

<^(o,...,o,A)( s )' one obtains: 

oo 

■km =(i+-r7(i^^ < i6) 



Its asymptotic behavior for z going to infinity is 

Vd (z) = l-^ z- 2 + 0(z-*) (17) 
whereas for z going to 0, one has 

%(*)~^|r^ (d>2), Vd (z)^-\nz (d = 2) 

These relations can be used for qualitative study of the partially reflected 
Brownian motion. For instance, one identifies the parameter A as a character- 
istic length scale of the problem: the magnitude of any distance (e.g., |s|) has to 
be compared with A. Interestingly, the asymptotic behavior (flTt for large z means 
that the function t\(s) is close to the harmonic measure density i^(o,...,o,a)( s )- 
Roughly speaking, for large |s|/A, the partially reflected Brownian motion started 
from the origin is qualitatively equivalent to the (simple) Brownian motion started 
from the point (0, 0, A). In other words, the partial reflections on the boundary 
lead to a spreading of the harmonic measure with characteristic scale A (see also 
relation STU ). The explicit analytical results can also be derived in presence of an 
absorbing barrier at a given height 15911 . 

The knowledge of the probability kernel Ta(s, s') brings an important infor- 
mation about the partially reflected Brownian motion in the (upper) half space. As 
an example, we calculate the probability T^(r) that the PRBM started from the 
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Figure 3: Land Surveyor Approximation: the total flux across the boundary can 
be approximately calculated when the mixed boundary condition dC/dn = C/A 
on a given irregular curve (on the left) is replaced by the Dirichlet condition C = 
on the coarse-grained boundary (on the right). The last one is obtained by replac- 
ing curvilinear intervals of length A by corresponding linear chords. 

origin is finally absorbed on the disk Bf~ x = { (xi, Xd) G K d : x\ + ... + 
x 2 d-i < r 2 , Xd = } of radius r centered at the origin: 

oo r/A 

r 2 r(-) r f x d ~ 2 dx 

p i( r) = J ds tA (.) = f ^ J te-'it J (18) 

B d-1 

This probability shows how far the partially reflected Brownian motion can go 
away after the first hit of the boundary. One sees that this function depends only 
on the ratio r/A, going to for small radii and to 1 for large radii. Again, the 
parameter A is the characteristic length scale of the problem. In two-dimensional 
case, V\ (A/2) is the probability that the PRBM is absorbed on the linear seg- 
ment of length A, centered at the origin (the first hitting point). The numerical 
calculation of the integral in (fl~8l gives Pa (A/2) ~ 0.4521, i.e., about half of the 
particles is absorbed on this region. In other words, the length of the character- 
istic absorption region (where half of the particles is absorbed) is approximately 
equal to A. It has been shown recently that this result is qualitatively valid for a 
large class of irregular boundaries 16011 . Roughly speaking, if the one-dimensional 
boundary (curve) has no deep pores (fjords) and its perimeter is large with respect 
to the scale A, then the curvilinear interval of length A, centered on the first hitting 
point, absorbs approximately half of the diffusing particles. 

This result can be considered as a first mathematical justification of the Land 
Surveyor Approximation (LSA) developed by Sapoval [19]. According to this 
approximation, a given one-dimensional interface (curve) can be coarse-grained 
with physical scale A in order to replace the mixed boundary condition [I — 
Ad/dn]C = by the Dirichlet condition C = (see Fig. [3l> . Some heuristic 
physical arguments allowed to state that the total flux across the irregular semi- 
permeable interface was approximately equal to the total flux across this coarse- 
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grained boundary with Dirichlet condition. This statement provided a simple but 
powerful tool to investigate Laplacian transport phenomena. The land surveyor 
approximation had been checked numerically IG3l Q], but not mathematically. 
The study of the partially reflected Brownian motion brings its justification and 
further understanding. Actually, the coarse-graining procedure generates the re- 
gions of length A, where about half of the particles is absorbed. The LSA is 
based on two simplifications which can be clearly explained in terms of diffusing 
particles: 

1. The Dirichlet boundary condition on the coarse-grained boundary means 
that all particles arrived to the characteristic absorption region are absorbed. 
This approximation does not take into account half of the particles which 
escaped this region. 

2. The linear chords, generated by coarse-graining, are deterministic regions. 
This approximation neglects the fact that the characteristic absorption re- 
gions should be centered at the random position of the first hit of the bound- 
ary. 

Although these simplifications seem to be rough, the numerical simulations show 
that the LSA reproduces the transport properties with good accuracy. However, 
this approximation has no any kind of small parameter which would allow to 
control its applicability. More accurate theoretical approaches will be discussed 
in Section |U 

One can go further by extending the Land Surveyor Approximation to the 
three-dimensional case, which still remains poorly understood. Indeed, the nu- 
merical calculation in 3D leads to Pa (A) ~ 0.461 1, i.e., about half of the particles 
is absorbed on the disk of radius A centered at the first hitting point. Consequently, 
if one finds a convenient cover of a semi-permeable irregular interface by disk-like 
sets of characteristic radius A, the LSA may be still valid, i.e., the total flux across 
a given interface would be approximated by the total flux across the perfectly 
absorbing coarse-grained interface (with Dirichlet condition). An accurate math- 
ematical formulation of this extension and its numerical verification present open 
interesting problems. 



3.3 Relation to the Dirichlet-to-Neumann Operator 

The construction of the partially reflected Brownian motion for a given domain VL 
requires the resolution of the stochastic differential equation ©, a quite difficult 
problem. Fortunately, many characteristics of this process, e.g., the spread har- 
monic measure density uj x> k(s), can be obtained in another way. This subsection 
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is devoted to the Dirichlet-to-Neumann operator and its relation to the partially 
reflected Brownian motion. 

Definition 3.7 For a given domain Q C M. d (d > 2) with smooth bounded bound- 
ary dQ, let u : Q U dQ — > R be a harmonic function with Dirichlet condition 
u = f, a function f being from the Sobolev space (in other words, u is 

the solution of the boundary value problem © )■ Applying the normal derivative 
to u, one obtains a new function g = du/dn belonging to the space L 2 (dQ) of 
measurable and square integr able functions. Then the operator ftA, acting from 
H 1 (dQ) to L 2 (dQ), which associates the new function g with a given f, is called 
Dirichlet-to-Neumann operator. 

It is known that the Dirichlet-to-Neumann operator M. is self-adjoint pseu- 
dodifferential operator of the first order, with discrete positive spectrum and 
smoot h eigenf unctions forming a complete basis in L 2 {dVl) Ipll I52L loll y&A l63i 
|^|65l|6a|67l|68|]. One can also define its resolvent operator T\ = [I + AAi]^ 1 , 
called spreading operator. This is an analytic operator function in the whole com- 
plex plane, except a denumerable set of points, C\{— /i^ 1 }. In particular, T\ is 
well defined for any positive A. 

Lemma 3.8 For any strictly positive A, the spreading operator Ta acts from 
L 2 (dQ) to L 2 (dQ) as a compact integral operator, 



where the kernel Ta(s, s') is given by ([77]). 

Proof . The probability kernel Ta(s, s') is a positive function satisfying the nor- 
malization: 



an 

since the partially reflected Brownian motion is conditioned to be finally absorbed 
on the boundary. Therefore, one obtains: 






an an 



where Stot is the total surface area of the boundary dQ. The integral operator Ta 
defined by the kernel T^(s, s') is a Hilbert-Schmidt operator and, consequently, a 
compact operator. 
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The boundary condition in lemma 13.41 can be written with the help of the 
Dirichlet-to-Neumann operator: 



that implies that the integral operator T\ coincides with the resolvent [I + AM.] 1 
of the Dirichlet-to-Neumann operator M. . □ 

This simple lemma creates a "bridge" between the partially reflected Brow- 
nian motion and the Dirichlet-to-Neumann operator. In particular, the relation 
(fTTT) for the spread harmonic measure density u} x \(s) can now be understood as 
application of the spreading operator T A to the harmonic measure density co x (s). 
Consequently, once the Dirichlet-to-Neumann operator is constructed for a given 
domain, one can calculate the density iv x> a(s) without solving the stochastic dif- 
ferential equations ©. 

The self-adjointness of the Dirichlet-to-Neumann operator allows one to apply 
efficient tools of the spectral theory. For example, one can rewrite the relation (ITTT) 
as spectral decomposition of the harmonic measure density on eigenfunctions V Q 
of the operator M : 



where ( ■ )l 2 denotes the scalar product in L 2 (dVL) space. The advantage of this 
relation is an explicit dependence on the physical parameter A. 

3.4 Examples 

In order to illustrate the underlying concepts, we consider several examples. 
3.4.1 Two-Dimensional Disk 

We are going to study the partially reflecting Brownian motion in a unit disk, Vt = 
{ x G M 2 : \x\ < 1 } (its boundary is a unit circle, dil = { x G M 2 : \x\ = 1 }). 

In this case, the harmonic measure density u> x (s) = cu(r, 6) is a function of 
two real variables: the distance < r < 1 between the starting point x G Vt and 
the origin, and the angle < 9 < 2n between directions onto points x and s G dVl 
from the origin. The harmonic measure density is known as Poisson kernel: 



[l + AM]T A (s,s') = 5(s-s') 





(19) 



u(r,8) 



27r(l-2rcos# + r 2 ) 



(20) 



Partially Reflected Brownian Motion. 



24 



The rotational invariance of the domain implies that the eigenbasis of the 
Dirichlet-to-Neumann operator M. is the Fourier basis, 

Add 

V a (9) = -= (aeZ) 
y In 

Taking Fourier harmonic as boundary condition, u(r = 1,9) = e iad , one finds 
a regular solution of the corresponding Dirichlet problem: u(r,9) = r^ a 'e ia0 . 
Since the normal derivative coincides with the radius derivative, one obtains the 
eigenvalues of the Dirichlet-to-Neumann operator: 

\i a = \a\ (a E Z) 

These eigenvalues are doubly degenerated (expect /i = 0). 

The spread harmonic measure density is given by relation JT9b : 



Wx,a(s) = u A (r,9) = — V j 
2ir ' 1 



j.\a\ giaO 



Aid 

(the scalar product of the harmonic measure density u x (s) and eigenfunctions 
V* (9) is shown to be equal to r' a l, with r = \x\). In the case A = 0, one retrieves 
the Poisson representation for the harmonic measure density (l20b just as required. 
The kernel of the resolvent operator T\ is 



9') 



■y 00 e ia(fl-8') 



2tt ^ 1 + Abl 



For the exterior problem, when Q = {x EM 2 : |x|>l}, one obtains exactly 
the same results. 



3.4.2 Three-Dimensional Ball 

The similar arguments can be applied for higher dimensions. For example, in the 
three-dimensional case, one considers the unit ball r2 = {sGlR 3 : \x\ < 1 }. 
The harmonic measure density is known to be 

1 -r 2 

u x (s) = u(r, 9) = — -jjz (s G dn) 

4n (1 - 2r cos9 + r 2 ) 1 

where r = \x\ < 1 is the distance between the starting point and the origin, 

and 9 is the angle between directions onto points x and s E dVl from the origin. 
This function can be expanded on the basis of spherical harmonics as 

oo I 

OJ. 
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The rotational symmetry of the problem implies that the eigenbasis of the 
Dirichlet-to-Neumann operator is formed by spherical harmonics Yi >m . A regular 
solution of the Dirichlet problem © in the unit ball can be written in spherical 
coordinates r, 9 and ip as 

oo I 

u(r, 9,<p) = Yl ^ m rl Y i,™{0, <p) 

1=0 m=-l 

where fi tm are coefficients of the expansion of a given boundary function / (Dirich- 
let condition) on the complete basis of spherical harmonics. Since the normal 
derivative coincides with the radius derivative, one obtains 

V on J an \ or J r=l i =0m =-l 
i.e., the eigenvalues of the Dirichlet-to-Neumann operator Ai are 

fH = l (I G{0,1,2,...}) (21) 

Note that the Z-th eigenvalue is degenerated n\ = (21 + 1) times. 

Interestingly 3 , the Dirichlet-to-Neumann operator Ai for the unit ball fl = 
{x G M 3 : \x\ < 1} coincides with an operator introduced by Dirac in quan- 
tum mechanics ll69ll . It is known that the hydrogen atom is described by three 
quantum numbers: the main quantum number n, the orbital quantum number I 
and magnetic quantum number m. Two last numbers are associated with indices 
of spherical harmonics. Thus, the Dirichlet-to-Neumann operator in the ball is 
apparently the orbital quantum number operator for the hydrogen atom. In par- 
ticular, the degeneracy of eigenvalues of this operator can be understood from the 
point of view of spin degeneracy. 

The spread harmonic measure density u> Xi \(s) and the spreading operator ker- 
nel T\(s, s') can be written explicitly as spectral decompositions on the eigenbasis 
of the Dirichlet-to-Neumann operator Ai as in the two-dimensional case. 

The eigenvalues of the Dirichlet-to-Neumann operator for rf-dimensional unit 
ball are still given by (l2TTi with degeneracy 

_ (2l + d-2) (Z + d-3)! 
ni ~ (d-2) (d-3)!Z! 

The exterior problem for fi = {a;GlR 3 : \x\ > 1} can be considered in the 
same manner. The harmonic measure density is u) x (s) = (1/r) u(l/r,6), with 

3 The author thanks Dr. S. Shadchin for valuable discussions on this relation. 
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r — \x\ > 1. Using the same expansion on spherical harmonics, one obtains 
Hi — I + 1 with I e {0, 1, 2, ...}. In particular, the lowest eigenvalue /i = 1 is 
strictly positive. This difference with respect to the spectrum fxi for the interior 
problem has a simple probabilistic origin: the Brownian motion in three dimen- 
sions is transient, i.e., there is a positive probability (equal to 1 — 1/r) to never 
return to the ball. Another explication follows from the theory of boundary value 
problems for elliptic differential operators: the exterior Neumann problem has a 
unique solution, while the solution of the interior Neumann problem is defined 
up to a constant. Consequently, the Dirichlet-to-Neumann M. operator should be 
invertible for the exterior problem that implies a simple condition for its eigen- 
values: ji a 7^ 0. On the contrary, M. is not invertible for the interior problem 
providing the condition that at least one eigenvalue is zero. 

4 Stochastic Approaches to Laplacian Transport Phe- 
nomena 

In this section, we return to the Laplacian transport phenomena, discussed at the 
beginning. First, we are going to introduce the notion of source, diffusing particles 
started from. The definition of the partially reflected Brownian motion will require 
only a minor modification. After that, a recently developed continuous approach 
will be presented with a special emphasis on its physical significance. Finally, 
we shall mention two other physical descriptions which can now be considered as 
useful approximations to the continuous approach. 

4.1 Notion of Source 

The description of the partially reflected Brownian motion given in the previous 
section does not involve a source, an important element for Laplacian transport 
phenomena. In this subsection, we are going to discuss the extension of previous 
definitions in order to introduce the source. As one will see, a minor modification 
will be sufficient. 

Throughout this subsection, we consider a bounded domain f2 with twice con- 
tinuous differentiable boundary composed of two disjoint parts, dfl and <9f2 , re- 
ferred to as working interface and source respectively 4 . In practice, the working 
interface and the source are well separated in space, therefore one may think about 
a circular ring as generic domain. 

Previously, the whole boundary of the domain had been considered as the working interface. 
For this reason, we preserve the same notation dVt for this object and hope that it will not lead to 
ambiguities. 
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As previously, one considers the reflected Brownian motion W t , started from 
any point i£flU dVt U 8VLq and reflected on the whole boundary dVt U dVt Q , the 
corresponding local time process £ t , and the stopping time defined in 13. 11 Let 
us introduce a new stopping time r as the first moment when the process W t hits 
the source dVt : 

t = mf{t > : Wt g dn } 

Then, the spread harmonic measure can be introduced for any subset A from Borel 
a-algebra B(dVt) (defined on the working interface alone!) as 

lu xA {A} = F{W n E A, T X A < r < oo} (22) 

We outline two distinctions with respect to the previous definition 13 .21 

• The measure is considered on Borel subsets of the working interface dfl 
only, whereas the reflected Brownian motion W t and the local time process 
£ t are defined on the whole boundary dVt U dVt Q . 

• There is a supplementary condition T\ < t providing that the process W t 
should be stopped (absorbed) on the working interface before hitting the 
source. 

Note that 1 — tu X: A{dQ} is the probability that the process W t started from a 
given point x E Vt hits the source dVt Q before its final absorption on the working 
interface dVt. 

One can easily extend the lemma l3~3l to this spread harmonic measure: 

Lemma 4.1 For any subset A from B(dQ) and fixed positive A, the spread har- 
monic measure uj x ^{A}, considered as a function of x, solves the boundary value 
problem: 



Au xA {A} = (xeQ), 



dn 



u xA {A} = I A (x) (xedQ) 

~u xA {A} =0 (xe dtt ) 

(23) 



Proof is similar to that of the lemma l3~3l The last condition holds since x e dQ 
implies r = 0. □ 

Corollary 4.2 Function C^x) = C (l — u x ^{dQ}) solves the boundary value 
problem (CJQ. 
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Proof is a direct verification. □ 

Consequently, a simple introduction of the source allows one to apply the 
previous description of the partially reflected Brownian motion to study Laplacian 
transport phenomena. Due to reversibility of the Brownian motion, one may think 
that (1 — (jj x h{d£i\)dx gives also the probability to find the partially reflected 
Brownian motion, started from the absorbing source, in dx vicinity of the point 
x E £1, under partially absorbing condition on the working interface <9f2. Note that 
such way of reasoning, being intuitive and useful, is quite formal. In particular, 
the (simple) Brownian motion started from the source returns to it infinitely many 
times with probability 1. If one really needs to define such a process, the starting 
point should be taken slightly above the source. 

Since the boundary d£l is supposed to be smooth, one can introduce the spread 
harmonic measure density uj x>a (s). In turn, the kernel of the spreading operator is 
defined as previously, T\(s, s') = uj SjA (s'). In particular, one retrieves the relation 
CD: 

Ux,a( s ) = / ds'uj xfi (s') T A (s',s) 



an 

where the harmonic measure density u x> o(s) is defined by relation d22l with A = 
0. 

The definition of the Dirichlet-to-Neumann operator can also be extended to 
domains with a source. For a given function / G H l (dVt), one solves the Dirichlet 
problem in the domain Vt: 

a n I U — f (X E dQ) 

Au = o (zefi), n ; ~ n \ 

K ' u = o [x e oll ) 

(in principle, one could consider another function on the source). For a given 
function / on <9f2, the Dirichlet-to-Neumann operator Ai, acting from ^(dtl) 
to L 2 {dVt), associates the new function g = du/dn on the working interface d£l. 
One can prove general properties of this operator and its relation to the partially 
reflected Brownian motion in a straight way. In particular, the spreading operator 
7a, defined by its kernel T^(s,s'), coincides with the resolvent operator [I + 
AAi]^ 1 . However, some normalization properties may be changed. In particular, 
for the probability kernel Ta(s, s'), one has 



J T A (s,s')ds' < 1 
an 



since the PRBM started from the working interface dQ can now be absorbed on 
the source. 
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4.2 Continuous Approach 

The stochastic treatment by means of the partially reflected Brownian motion 
brings the solution to the problem ([l]-0) describing Laplacian transport phenom- 
ena: C A (x) = Co(l — ui x>A {d£l}) (see corollary 14.21) . One can go further using 
the close relation to the Dirichlet-to-Neumann operator 1E9I kdl . According to 
the lemma |4~T1 the density u> XtA {dQ}, considered as a function of x, solves the 
boundary value problem: 



Au xA {dtt} = (x e fi), 



dn 



tu xA {dtt} = (x E dn ) 



The restriction of the function uj x A {d£l} on d£l can be written with the help of 
the Dirichlet-to-Neumann operator M. as 

LO SjA {dQ} = [(I + AM)- 1 !]^) = [T A l){s) 

where 1 stands for a constant (unit) function on the working interface. 

One defines then the flux density 4>\(s) across the working interface OVl: 

Ms) = D —(s) = -DC ^ {s) 

Since the normal derivative of a harmonic function can be represented as the ap- 
plication of the Dirichlet-to-Neumann operator to the restriction of this function 
on the boundary, one writes 

A (s) = DC [Mu StA {dQ}](s) = DC [MT A l](s) 

(the sign is changed due to particular orientation of the normal derivative). Taking 
A = 0, one finds <p (s) = DC [Ml)(s) and finally 

0aO) = [T A <j) Q ](s) 

The transport properties of the working interface can be characterized by a 
physical quantity called spectroscopic impedance. We remind that the impedance 
of an electric scheme is defined as the tension applied between two external poles, 
divided by the total electric current passing through. The formal analogy between 
the electric problem and the diffusive transport, discussed in Section d leads to 
a natural definition of the impedance in our case as the concentration C on the 
source d£l divided by the total flux across the working interface <9fh 

Z ce ii(A) = f , ° : 
J ds(p A {s) 

an 
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Taking A = 0, one deals with a purely absorbing interface d£l: any particle 
arrived to d£l is immediately absorbed (without reflections). In other words, such 
interface has no resistance for passage across it. Consequently, the impedance 
■^ceZz(O) represents the "access resistance" by the bulk: the possibility that the 
Brownian motion can return to the source without hitting the working interface. 
The resistance of the working interface alone can thus be characterized by the 
difference between Z ce u(A) and Z ceH (0), called spectroscopic impedance: 

Z sp (A) = Z cell {A) - Z cell (0) 

Using the simple identity 

G)((0o-<M-1) L2 = ^ (0a-0o) l2 (24) 
one writes the spectroscopic impedance as 



Z sp {A) 



A (<Pa-<Po) L 2 



D (0 A -1) L2 (0 O -1) L2 
Applying again the identity (041) . one finds a more convenient form: 

1 



Z sp (A) 



1 



Z(A) Z cell (0) 

The new function 

^(A) = ^(T A ^-^) L2 
can be called effective impedance, where 

j.ht _s 



•1) 



is the normalized flux density towards the perfectly absorbing working interface 
<9f2. Finally, the spectral decomposition of the spreading operator Ta on the ba- 
sis of eigenfunctions V a of the Dirichlet-to-Neumann operator M. leads to the 
important relation for the effective impedance: 

Z ^ = £ E rri^r Fa = W • ■ v :), 2 ( 25 ) 

a 

This relation presents the central result of the continuous approach developed 
in lEoll . Let us briefly discuss its physical meaning. The spectroscopic impedance 
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Z sp (A) or, equivalently, the effective impedance Z(A), is a physical quantity that 
characterizes the transport properties of the whole working interface. More impor- 
tantly, this quantity can be measured directly in experiment (e.g., in electrochem- 
istry). On the other hand, the local transport properties of the working interface 
are described by the single physical parameter A, being related to the membrane 
permeability W, the electrode resistance r or the catalyst reactivity K (see Sec- 
tion [T]). Varying the parameter A, one changes the local transport properties at 
each boundary point and, consequently, the whole linear response of the working 
interface. At first sight, one may think that an increase of the local boundary re- 
sistance would imply a proportional increase of the whole boundary resistance, 
i.e., Z(A) ~ A. This reasoning, being true for a planar surface, becomes invalid 
in a general case due to geometrical irregularities and related screening effects. In 
fact, an irregular geom etry modifies considerably the linear response of the work- 
ing interface lElTH SH H H EJ B. The boundary value problem 



(U]-|3l), describing Laplacian transport phenomena on the average, allows formally 
to study such geometrical influence. Practically, however, this is a very difficult 
problem. In contrast, the continuous approach provides an efficient tool to carry 
out these studies both in theoretical and numerical ways. In particular, the relation 
(l25l) makes explicit the impedance dependence on the local transport properties 
(parameter A) and allows one to identify contributions due to the physics and due 
to the geometry, originally involved in the problem in a complex manner. In other 
words, whatever the physical problem (diffusion across semi-permeable mem- 
branes, heterogeneous catalysis or electric transport), the geometry enters only 
through the spectral characteristics of the Dirichlet-to-Neumann operator Ai: its 
eigenvalues p, a and the spectral components F a of the normalized flux density 
0q(s) on the basis of its eigenfunctions V a (s). 

The other important meaning of the relation (l25t can be outlined if one con- 
siders the inverse problem IG9I1 : what is the most available information that one 
can retrieve from a measurement of the spectroscopic impedance of an unknown 
working interface? The mathematical response can be given immediately if one 
rewrites (l25t as Laplace transform of the new function C(A): 



oo 

£(A) = i J dAe- A / A C(A), C(A)e=2>« 



Under assumption to be able to measure the impedance with an absolute preci- 
sion, one can reconstruct the function £(A) and, consequently, the set of charac- 
teristics F a } which may thus be called harmonic geometrical spectrum of 
the working interface. The hierarchical structure of this spectrum for self-similar 
boundaries has been recently investigated iFoll . Many interesting properties of the 
function £(A) remain poorly understood. 
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In the two following subsections, we are going to discuss some aspects of the 
semi-continuous and the discrete descriptions of Laplacian transport phenomena. 
These approaches are based on more intuitive notion of partial reflections on the 
boundary. Since these descriptions turn out to be approximations to the continuous 
approach, we do not present the circumstantial details. 



4.3 Semi-continuous Approach 

Halsey and Leibig gave the first theoretical description of the electrolytic double 
layer response with emphasis on electrochemical applications W29I1 . This descrip- 
tion, involving the Green function of the electrolytic cell, can be reformulated in 
the following stochastic language (for details, see |39]). For a given domain 
with smooth bounded boundary dVl, one considers the Brownian motion started 
from a point iGfl. When the diffusing particle hits the boundary at some point 
s, two complementary events may happen: 

• with probability e, the Brownian motion is reflected to the interior bulk point 
s + an(s), slightly above the boundary (here n(s) is the unit normal vector 
to the boundary at point s, a is a small positive parameter); the Brownian 
motion continues from this point; 

• or, with probability 1 — e, the Brownian motion is terminated at this point s 
(absorbed on the boundary). 

This stochastic process is continued until the absorption on the boundary and can 
be called Brownian motion reflected with jump. Two new parameters, the jump 
distance a and the reflection probability e, are related to the given physical length 
A Hi: 

£ = TTT7T\ (26) 
1 + (a/A) 

Now, one can calculate the probability uj^ A (s)ds that this process is finally 
absorbed in ds vicinity of the boundary point s. Since the motions before and 
after each reflection are independent, this probability can be obtained as the sum 
of probabilities to be absorbed after 0, 1, 2, ... reflections: 



£h(s)ds = [u x (s)ds] (l-e) + J [w<c(si)cisi] e [uj sl+an ( Sl )(s)ds] (1 - e)+ 

+ J J [u x (si)dSi] £ [uJ Sl+ an{s 1 ){s2)ds 2 ] E [u S2+an{s2 ) (s)ds] (1 - e) + 



an an 



For example, the third term represents the probability to hit the boundary in ds^ 
vicinity of the point s 1: to be reflected to the neighboring point si + em(si), to hit 
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again the boundary in ds% vicinity of the point s 2 , to be reflected to the neighboring 
point s 2 + an(s2), to hit the boundary for the last time in ds vicinity of the point 
s, and to be finally absorbed. Introducing the integral operator Q^ a \ acting from 

L 2 (dn) toL 2 (dn) as 

[Q (a) f](s) = J ds'f(s') co s , +an(sf) (s) 

one rewrites the previous sum as the application of the new integral operator 
to the harmonic measure density lo x (s): 

oo 

= [tP^](s) with if = (1 - e) J2( £ Q (a) ) k ( 27 ) 

fc=0 

What happens when the jump distance a goes to 0? Hitting the boundary, the 
Brownian motion will be reflected to interior points lying closer and closer to the 
boundary, i.e., displacements of the Brownian motion between two serial hits are 
getting smaller and smaller. At the same time, the reflection probability e tends 
to 1 according to relation (l26l) . i.e., the average number of reflections increases. 
Indeed, the distribution of the random number J\f of reflections until the final 
absorption is simply 

F{Af= n} = (1 -e) e n (28) 

implying that the average number E{7V} = e(l — e)^ 1 goes to infinity. Does a 
limiting process exist? The situation is complicated by the local choice between 
reflection and absorption: at each hitting point, the motion can be absorbed with 
vanishing probability 1 — e. In order to overcome this difficulty, one can consider 
this process from a slightly different point of view. Actually, one can replace the 
local condition of the absorption (with probability 1 — e) by its global analog: the 
process is absorbed on the boundary when the number of reflections exceeds a 
random variable J\f distributed according to the geometrical law d28t . Evidently, 
this modification does not change at all the properties of the process. At the same 
time, we gain that the condition of the absorption becomes independent of the 
Brownian motion between serial hits. As a consequence, one can consider the 
corresponding limits (as a — > 0) separately. So, the Brownian motion reflected 
with jump should tend to the reflected Brownian motion as the jump distance 
a vanishes. This motion, however, is conditioned to stop when the number of 
reflections exceeds the random variable J\f. Since the average number E{A/"} goes 
to infinity in the limit a —> 0, it is convenient to consider a normalized variable 
X = aM obeying the following distribution: 

oo 

P{ x >\} = F{Af > A/a} = F {N = n}~ s [x/a] ~ exp[-A/A] (29) 

[A/a] 
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(the last equality is written with the help of (l26t for a going to 0). Since the 
number of reflections on jump distance a, multiplied by a, tends to the local time 
process according to Levy's formula ©, the previous condition of absorption can 
be reformulated: the motion is absorbed when its local time process exceeds a 
random variable distributed according to the exponential law (l29l) . One thus con- 
cludes that the Brownian motion reflected with jump should tend to the partially 
reflected Brownian motion defined in 13.11 

The above analysis, presented as a sketch (without proofs), does not pretend 
to a mathematical rigour. It may be considered rather as a possible justification 
which can be brought for the semi-continuous approach if necessary. In particular, 
one can demonstrate that the density ^ a A (s), given by relation (l27l) . tends to the 
spread harmonic measure density as the jump distance a goes to 0: 

u xA (s) = limu^s) 

This relation may be useful for numerical computations (in particular, it was ap- 
plied in [60]). Similarly, the integral operator should converge to the spread- 
ing operator T A as a — > 0. Calculating the geometrical series in (|27l) and repre- 
senting as 

Ti a) = (l-e)(l-eQ^Y 1 =(l + A 1 -^^) 1 

one obtains the following approximation for the Dirichlet-to-Neumann operator: 

/ - Q (a) 
M = lim 

a^Q a 

Again, this relation may be useful for the numerical computation of this operator. 

The advantages of the semi-continuous approach are based on an apparent 
intuitive meaning of partial reflections on the boundary. Moreover, this approach 
provides even a more realistic description of physico-chemical processes at micro- 
scopic level. For example, if one considers the partially reflected Brownian motion 
started from a boundary point, the number of hits of the boundary is infinite for 
any moment t > that sounds impossible for real physical species. The keypoint 
is that, for diffusion across a semi-permeable membrane or heterogeneous reac- 
tion on a catalytic surface, the description by the boundary value problem dU-Hl) 
cannot be justified on length scales less than the mean free path of diffusing parti- 
cles. Since the continuous limit a — > requires such non-physical scales, it is not 
surprising that the limiting process (the PRBM) presents some irrealistic proper- 
ties from the physical point of view. The similar limitation happens for the electric 
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transport problem for which the smallest physical scale is given by the thickness 
of the double layer, being close to the Debye-Hiickel length [26,|22|]. Evidently, 
this remark does not devaluate the efficiency of the continuous approach based 
on the partially reflected Brownian motion. On the contrary, the mathematical 
rigour of this approach justifies the semi-continuous description and simplifies its 
study by introducing the Dirichlet-to-Neumann operator. However, when dealing 
with a mathematical description of a physical problem, one should take care that 
deduced consequences do not go beyond the ranges of the model. 

The capabilities of the semi-continuous approach are essentially limited by the 
fact that the governing operator is not self-adjoint (the function u) s+an ( s \(s') 
is not symmetric with respect to the permutation of s and s' except specific cases). 
As a consequence, one cannot develop the spectral decomposition (l25l) of the 
impedance. In particular, there is no possibility to distinguish contributions from 
different eigenmodes. Although the operator is defined naturally by the har- 
monic measure density, it does not provide a proper description of the problem as 
it was done with the Dirichlet-to-Neumann operator. 

4.4 Discrete Approach 

Another stochastic approach to Laplacian transport phenomena was developed 
by Filoche and Sapoval [2 1 J . The main idea is to model the partially reflected 
Brownian motion by lattice random walks with partial reflections on the bound- 
ary. Actually, one discretizes a given domain f2 by <i-dimensional hypercubic 
lattice of mesh a and considers the following stochastic process: started from a 
remote source, a random walker jumps to a neighboring site at each step with 
probability (2d) -1 . When the walker arrives to a boundary site, it can be reflected 
to its neighboring site (belonging to the bulk) with probability e (and the motion 
continues), or it can be absorbed with probability (1 — e). The motion continues 
until the final absorption on the boundary, or the return to the source. One can 
show 12211 that the discrete parameters a and e are related by the expression (l26l) 
involving the continuous physical parameter A. 

In the discrete description, the harmonic measure density is replaced by the 
distribution of hitting probabilities (Po)j on boundary sites j, (simple) random 
walks being started from a remote source. Let Qj k denote the probability to 
arrive to the boundary site k starting from the boundary sites by a random walk 
in the bulk without hitting the boundary or the source during the walk 5 . One can 
thus calculate the distribution of probabilities (Pa) j to be finally absorbed on the 

5 We use the same notation for the integral operator in semi-continuous approach and 
for the matrix of these probabilities since they have the same meaning and even may be used to 
approximate each other. 
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boundary sites j, when random walks with partial reflections are started from a 
remote source. Indeed, the Markov property of this process allows to calculate 
(Pa)j as me sum °f contributions provided by random trajectories with 0, 1, 2, ... 
reflections before the final absorption: 



0)j{ 



^)+E( p o)^^£-( 1 - £ )+EE( p o)^^i a i 2 ^Si( 1 - £ )+ 



fel k2 



(we remind that e = (1 + a/A) -1 ). For example, the second term represents the 
product of the following probabilities: to hit a boundary site k\, to be reflected to 
its neighboring site, to arrive to the boundary site j, and to be finally absorbed on 
it. If one considers Po and Pa as vectors and QW as matrix, the summation over 
intermediate sites fci, fe, ... can be understood as matrix product: 



;i- £ )EK> (a) ) n 



71=0 

i.e., the distribution of absorption probabilities (Pa),? is obtained as the applica- 
tion of a linear operator, depending on and A (or e), to the distribution of 
hitting probabilities (Po)j- The symmetric matrix Q represents a self-adjoint 
operator, called Brownian self-transport operator. Using the normalization prop- 
erty \Q^ a > | < 1 and relation (1261) between A and e, one obtains: 



71(a) -p 



i(a) 



J + A 



I-Q 



(a) 



The operator , depending on the lattice parameter a, is called ( discrete ) spread- 
ing operator. The previous relation, written explicitly as 



B P o)47l a) ) 



allows one to separate random trajectories in two independent parts: 

• the random walker started from a remote source arrives to the boundary site 
k (first factor); 

• it continues the motion with partial reflections until the final absorption on 
the boundary site j (second factor). 

One concludes that the absorption probabilities (Pa) j provide a discrete ana- 
log of the spread harmonic measure density, while the matrix (T^) is a dis- 
crete analog of the kernel Ta(s, s') of the spreading operator T\ = [I + AA4]~ L . 
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In particular, the bounded operators (I — Q )/° can be understood as discrete 
approximations of the Dirichlet-to-Neumann operator M. (in resolvent sense). As 
for the semi-continuous approach, we do not furnish the corresponding proofs (see 
ll39ll for more details). 

The advantage of the discrete description with respect to the semi-continuous 
approach is based on the fact that the Brownian self-transport operator Q^ a > and, 
consequently, the (discrete) spreading operator are self-adjoint. This property 
allows to employ all the machinery of the spectral theory in order to express the 
physical characteristics of Laplacian transport through eigenmodes of this opera- 
tor in an explicit way. For example, the spectral decomposition (l25t can be writ- 
ten in the discrete case. Such decompositions have been used to study Laplacian 



transport towards irregular geometries [22, 39]. Moreover, the discrete descrip 



tion suggests at least two different ways to study the problem numerically: direct 
Monte Carlo simulations and discrete boundary elements method. 

The discrete description, being intuitively the most simple and useful, may 
lead to mathematical difficulties when one tries to proceed the continuous limit 
a going to 0. Although the partially reflected Brownian motion is the natural 
limit of random walks with partial reflections, its rigorous demonstration, in our 
knowledge, is not yet realized in details. The interested reader can find more 
information on this topic in iFl . 72 . 73 . 74 . 75 . 



5 Conclusion 

The application of stochastic processes to represent the solution of boundary value 
problems is well known and wide used. In particular, Monte Carlo simulations are 
generally based on this concept. In this paper, we gave a brief overview of Lapla- 
cian transport phenomena in different scientific domains (e.g., physics, electro- 
chemistry, chemistry, physiology) and related stochastic approaches to describe 
them. The most attention has been paid to the recently developed continuous ap- 
proach based on the partially reflected Brownian motion. This stochastic process 
can be thought as rigorous mathematical description for random trajectories of 
diffusing particles hitting a semi-permeable interface, in comparison with more 
intuitive physical descriptions by semi-continuous and discrete approaches. The 
partially reflected Brownian motion turns out to be the natural limit of the Brow- 
nian motion reflected with jump (semi-continuous approach) and of the lattice 
random walks with partial reflections (discrete approach). 

The profound relation between the partially reflected Brownian motion and 
the spectral properties of the Dirichlet-to-Neumann operator Ai are shown to be 
useful for practical purposes. In particular, the kernel of the resolvent operator 
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Ta — [I + AM]' 1 gives the probability density T A (s, s') allowing to reconstruct 
the spread harmonic measure uj Xj a- Moreover, the spectral decomposition on the 
complete basis of the Dirichlet-to-Neumann operator eigenfunctions leads to the 
explicit analytical formula for its density. Consequently, the use of the operator 
M. is an efficient way to study different probability distributions related to the 
partially reflected Brownian motion. 

The spectral decomposition of the spectroscopic impedance, characterizing 
the linear response of the whole working interface, leads to an explicit analytical 
dependence on the physical parameter A allowing to identify physical and geo- 
metrical contributions which were involved in a complex manner. The harmonic 
geometrical spectrum of the working interface contains the complete information 
about its transport properties. The combined use of stochastic characteristics of 
the partially reflected Brownian motion and spectral properties of the Dirichlet- 
to-Neumann operator opens encouraging possibilities for further understanding 
various physical and chemical transport processes in nature. In this light, a more 
profound mathematical analysis of these objects seems to be an important per- 
spective for the present study. 
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